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Abstract 

"S Stark-chirped rapid adiabatic passage (SCRAP) is an important technique used for coherent quantum 

^ controls. In this paper we investigate how the practically-existing dissipation of the system influences on 

'— ' the efficiency of the passage, and thus the fidelities of the SCRAP-based quantum gates. With flux-biased 

►^ Josephson qubits as a specifical example, our results show clearly that the efficiency of the logic gates 

QQ implemented by SCRAP are robust against the weak dissipation. The influence due to the non-adiabtic 

cn 

■^ transitions between the adiabatic passages is comparatively significantly small. Therefore, the SCRAP- 



based logic gates should be feasible for the realistic physical systems with noises. 
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Over recent years, quantum computation has attracted much attention partly because the discov- 
ery of quantum algorithm for specific problems provides a tremendous speedup in computation, 
compared to a classical computer [HI El- A critical prerequisite for building a quantum computer 
is to perform the basic single- and two-qubit gates with high fidelity above certain threshold lev- 
els m^- 

A typical ingredient in quantum computing is the coherent transfers of the population between 
the qubit states. Basically, there are two approaches to realize the population transfers between 
two selected quantum states; one makes use of the Rabi oscillations and the other is based on pop- 
ulation passages [|51[6l|. For Rabi oscillations, the transfer efficiency is strongly dependent on the 
precisely-designed duration of the applied pulse. On the other hand, the logic gates implemented 
via population passage techniques, such as the stimulated Raman adiabatic passage (STIRAP) [|7]| 
and the Stark-chirped rapid adiabatic passage (SCRAP) [8J, are evolution-time insensitive and 
thus robust against the imperfections of durations of the applied pulses. Until now, most of pop- 
ulation passage schemes to implement the quantum computation are based on the pure quantum 
systems, but their practically-existing dissipative effects (e.g., spontaneous emissions, phase re- 
laxations and the outsides from the system) have not be exactly treated. The fidelities of the logic 
gates for the realistic quantum computing demonstrations are particularly important, therefore, it 
is necessary to investigate how the practically-existing dissipation influences on the efficiencies 
of the population passages and consequently the fidelities of the relevant logic gates for quantum 
computing. 

Usually, the dynamics of an open quantum system can be described by two approaches AH: the 
master equation for the reduced density matrix and the Heisenberg-Langevin equation by intro- 
ducing environment noise operators. Additionally, a relatively-simple approach, i.e., dissipative 
Schrodinger equation with a non-Hermitian Hamiltonian describing the damping, is also utilized. 
In this approach the environment effects are considered simply by phenomenologically introduc- 
ing certain non-Hermitian terms in the Hamiltonian of the system. Then, the dynamics of the 
treated open system can still be described by the usual Schrodinger equation. Indeed, this idea 
has been utilized to investigate the dissipative effects in the well-known STIRAP for three-state A 
atomic systems [fTOl . wherein the decay rate T of the intermediate state |2) is served as the main 
source of the dissipation during the population transfers from the state |1) to the target state |3). 
The damping of the transfer efficiency with T shows different behaviors, i.e., exponential at small 
T and polynomial at large T. This feature provides a realistic STIRAP scheme for three-state A 



atomic systems in the presence of decay of the intermediate state. For the present two-state sys- 
tem, the decay of the excited state is the dominant dissipative source which mainly decreases the 
transfer efficiency of the SCRAP [jTTlfTll . In this brief report, we examine how this dissipation 
influences the fidelity of the SCRAP-based logic gates. For the simplicity, we treat the dissipation 
related to the excited state of the qubit by adding an imaginary part to the relevant diagonal term 
of the non-dissipative Hamiltonian. Our proposal is demonstrated specifically with the flux-biased 
Josephson qubits, but can also be applied to the other driven solid-state qubit systems. 

Without loss of the generality, we assume that the two-state system is well prepared initially, at 
time to = — C)0, in the ground state |0). Our end is to maximize the final population Pi(oo) of the 
target state 1 1) in the presence of the decay of 1 1) . Originally, without dissipation of the two- state 
system the desired transfer can be precisely implemented by means of the standard SCRAP fT3L 
151 . However, due to various practically-existing noises, such a transfer should be influenced. As 
a simplified model the state decay of a quantum system can be generically described by adding 
a loss rate T (> 0) to its non-lossy Hamiltonian, as a negative imaginary part to the correspond- 
ing diagonal term. As a consequence, the time evolution of the probability amplitudes for the 
dissipative driven two-level system can be expressed by the equation llTOl[TT]| 

Co\hf Q{t) \ fCo\ ^^^ 

Ci j 2 \^ ^(^) 2A(t) -2tT J \Ci j ' 

Here, fi(t) is the Rabi frequency coupling the levels of the two-state system, A(t) is relative to the 
pulse chirping the excited level and Co(t) and Ci(t) are the probability amplitudes related to the 
states |0) and |1), respectively. 

To analyze the progress of SCRAP in the presence of state decay, we define the adiabatic states 
1+) = sin6'(t)|0) + cos6'(t)|l) and |— ) = cos^(t)|0) — sin6'(t)|l), which are the instantaneous 
eigenstates of the Hamiltonian in Eq. ([T]) with F = 0. Here, the mixing angle 9{t) is defined as 
9{t) = arctan[r2(t)/A(t)]/2. In the basis defined by the adiabatic vectors |+) and |— ), equation 
([T]) can be written as 

a+\ 1 I e+{t) 2i9 
a j 2 1 _2^0 £_(t) 






— iT cos^ 9 iT sin 9 cos 9 

+ 1 . . 2. I I ' h (2) 

iV sin 9 cos 9 —iV sin 9 



with £±(t)=A(t) ± ^yA'^{t) + n^{t). 



Obviously, the off-diagonal elements in Eq. (|2]) result in the coupling between two adiabatic 
states 1+) and |— ) (i.e., the passage paths for the desired population transfers). For the ideal case 
without state decay, i.e.,r = 0, the desired adiabatic transfer can be implemented by properly 
designing the applied pulses to satisfy the condition: 9 = 0, i.e., the adiabatic condition llT6l[T7ll 

^ m)dA{t)/dt-A{t)dnit)/dt\ 

However, the second term in Eq. ([2]) shows that the damping of SCRAP corresponds to two ways, 
one is the decay of the adiabatic passage paths described by the diagonal elements and the other 
is the transition damping described by the nonzero off-diagonal elements. For a counterintuitive 
pulse sequence with initial state |0) (at time t = — oo, ^ = 0, then 9 = n/2 at t = oo), the transfer 
progress (which transfers the population from state |0) to state |1)) goes along the adiabatic path 
|— ) with the decay rate T sin^ 9. Along this adiabatic passage, the final population of state 1 1) is 

Pf ^ exp ( -2r f sin^ 9{t)dt\ , (4) 

while if the system is initially prepared at the state 1 1), then the population is transferred along the 
adiabatic passage |+) (with the decay rate F cos^ 9) to the ground state |0). The final population 
of the state |0) reads 

P^' ^ exp ( -2r I cos^ 9{t)dt\ . (5) 

Absolutely, the non-adiabatic transition between the states |— ) and |+) may also lead to the losses 
of Pf^ and Pq\ As the dissipation is irreversible, the population transfer may be significantly 
destroyed by the strong dissipation. 

In what follows, we investigate specifically how the dissipation of system influences the fidelity 
of the SCRAP-based quantum logic gates. Our discussion is based on the SCRAPs in flux-biased 
Josephson qubits, but can be easily generalized to other physical systems. For operational sim- 
plicity, here linear Stark pulses, rather than the previous gaussian Stark pulses [fT6l . are applied to 
the qubits. 

The quantum behavior of a flux-biased Josephson junction has been described in detail else- 
where [[T8l[T9ll . The Hamiltonian of the system is 

^ ^ 2m V 2A 



Here, the pump pulse lac = ^{i) cos{ujiot) is used to couple the qubit states and the Stark pulse Idc 
is applied to chirp the qubit's transition frequency Uiq. Also, $o = h/2e is the flux quantum, Ej = 
/o$o/2vr is the Josephson energy, and A = 27r/oL/$o, m = Cj[<l>o/(27r)]^, 0^ = 2ttI^qM/^o. 
Consequently, the Hamiltonian of the driven qubit (with decay rate F) in the interaction picture 
can be expressed as 







-t^^oi 



where 5i, = ^\5\j), i,j = 0, 1, k = ^{t)/2 and Ai(t) = MJrf,(t)(5n - 5oo)/L. 



I^kS^o -gAi(t)-znr 



(7) 
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FIG. 1: (Color online)The population transfer without dissipation for implementing the single-qubit gate 
with a flux-biased Josephson junction. In (a) the two pulses are designed with a linear forms ldc{^) = O.lt A 
and ,^(t) = —1.88 nA (—3.5 ns < t < 3.5 ns, else where ^(t) = V/m). With this pulse sequence, the 
system initially prepared in the state |1) completely passages to the state |0). The corresponding adiabatic 
parameter is shown in (b). 

When r = 0, i.e., for the ideal system without dissipation, we show in Fig. 1(a) that the 
single-qubit gate, i.e., the qubit inversions, can be realized by using a linear pump pulse /„c and 
a Stark pulse I^c to implement the desirable population transfer between the qubit states. It is 
shown that, under the counterintuitive pulse sequence (the applied Stark pulse I^c precedes the 
pump pulse but turns off first), the qubit inversion is realized along the adiabatic passage | — ) (with 
100% probability). Fig. 1(b) exhibits that the adiabatic parameter ri is fairly smaller than 1. This 
implies that the above progress for population transfers is really confined in the adiabatic region. 
Note that the desired population inversions are finished within a relatively-short time interval, i.e., 



Ti = 20 ns, which is really rapid compared to the typical decoherence time (e.g., 0.3 /is 11201 '). 

Now, let us consider how the dissipation of the system influences the above qubit inversions. 
For the convenience, we introduce a dimensionless decay rate 7 = TT [fTOl , where T is the 
characteristic width of the driving pulses. Then, the dissipation of the system can be divided 
into three regions; (i) weak dissipation (7 <^ 1), (ii) strong dissipation (7 ~ 1) and (iii) very 
strong dissipation (7 ^ 1). In Fig. 2 we show how the population probability of the target state 
varies with the decay rate 7 and the evolution time t for the applied counterintuitive sequence 
pulses. Specifically, Figs. 2(a) and 2(b) illustrate the population passage from the initial state |0) 
to the target state |1) along the adiabatic passage | — ); while Figs. 2(c) and 2(d) are relative to 
the population transfer from the state |1) to the state |0) along the adiabatic passage |+). The 
time-dependent population probabilities of the target state are calculated by Eqs. (|4]), Q. As 
a comparison, we also provide the relevant results by directly solving the Schrodinger equation 
with Hamiltonian Q. Here, we assume the qubit is in the initial state at time to = — 10 ns, the 
passage transfer is finished at time t/ = 10 ns, and the system is in the superposition state during 
the time tb = —3.5 ns to tm = 3.5 ns. It is shown that the adiabatic approximation made for 
delivering Eqs. (|4]) and ^ works well. The above numerical results show clearly that: (i) in the 
weak dissipation region, i.e., 7^1, the dissipation can be really neglected, and the efficiency of 
the population transfer is sufficiently high (almost 100%); in the strong dissipation, i.e., 7 ~ 1, 
the population efficiency is lower than 1; the final population may completely destroyed with a 
very strong dissipation 7^1. (ii) For the population passage from the state |0) to the state |1), 
the transfer probability decreases as an exponential function exp[— 2r(t — t„i)] after the passage 
region t > tm', while, for the passage from the state |1) to the state |0), the loss of the population 
can still be described by an exponential function ex'p[—2T{t — tb)] before the passage region t < tb- 
(iii)The non-adiabatic transition weakly influences the population transfer, and the dissipation of 
the system is mainly from the decay of the adiabatic passage paths. 

To be more thorough, we investigate how the dissipation influences the SCRAP-based two- 
qubit gate with two capacitively-coupled flux-biased Josephson qubits. For the simphcity, here the 
two junctions are assumed to be identical and thus two qubits possess the same energy structure. 
Originally, the two-qubit gate can be implemented also by the adiabatic population passages [|8l [T6ll 



(2) 
via applying a controllable dc current I^^' to chirp the second qubit. Considering the practically- 
existing decay of the excited state of the qubits (with the same decay rate T for simplicity), the 
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FIG. 2: (Color online) Population transfer with different decay rate for single-qubit gate. The pulses used 
to implement the adiabatic passage are the same linear pulses with the counterintuitive sequence for single- 
qubit gate discussed without dissipation. The population probability with initial state |0) at time to = 
— 10 ns varies with 7 = TT (T = 2 x 10^^) described by (a), while (c) is relative to the initial state |1). 
The red lines both in (a) and (c) are obtained by numerical solution of the Shrodinger equation related to 
Eq. 0, and the dashed green line in (a) and (c) is the analytical results from Eq. Q and Eq. (|5) respectively. 
Finally, (b) and (d) give the dynamics of the population for the different 7 and the evolution time t with 
initial state |0) and |1), respectively. 



Hamiltonian of such a driven two-qubit system can be simply expressed as 

/Aoo \ 

Aoi - ihr ^01 

fiio Alo - ihr 

\ All - 2ihr ) 



Hiit) 



(8) 



with 



Aoo - - -^—r^dc l*J^oo + [-T-) T^PooPoo > 
ZnL <Po Cm 

_ ^'^"^0.(2).^. , r^TT 2 1 (1) (2) 

27rL <Po Cm 

All - - ^;^/.c (i)^ii + (^) ^PiiPii ' 



•'m 



and 



o -o - r^^^2 1 m f2) 
•^0 C„ 



•'m 



where Cm = Cji^+C)/C (C is the effective coupling coefficient) represents the interaction between 
two qubits and pij = -ih{i\-§^ \j), plj = -pi^ . 

Still, one can easily check that the populations of |00) and |11) of the present two-qubit system 
are always unchanged, and the population transfer only occurs between the states |01) and |10). 
So the dynamics of the two qubits can be limited to a 2 x 2 subspace generated by the states |01) 
and 1 10). In absence of the dissipation, i.e., F = 0, Fig. 3(a) shows that the population transfer can 
be easily achieved between the states |01) and |10). Fig. 3(b) displays that the maximum value of 
the adiabatic parameter r/ during such a passage is about 0.14. Thus, the usual i-SWAP gate has 
been realized by the adiabatic SCRAP technique. 

In Figs. 3(c) we investigate how the dissipation influences the population transfer from the state 
|01) to 1 10) for a defined passage time interval T2 = 400 ns. It is shown that results by numerically 
solving the Schrodinger equation with for the Hamiltonian ([8]), and those by analytically solving 
the evolution within the subspace with the reduced Hamiltonian (defined by the adiabatic vectors 
|+)and|-)) 

HKt) ={ , (9) 

\^ e^-ihV j 

with e-t = (Aio — Aqi ± A/4n^7+TAio^^^oi)^)/2 are consistent. Obviously, the dissipation 
of the two-qubit operation is not relative to the non-adiabatic transition between the two passage 
paths 1+) and | — ). In Fig. 3(d) we depict how the transfer probability depends on the dissipation 
parameter 7 and the evolution time t. We can see from the figure that, for the sufficiently- weak 
dissipation (typically for 7 < 0.1) the passage time could be set as a sufficiently-long interval, 
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FIG. 3: (Color online) Population transfer for the two-qubit gate with a Stark pulse I^^ = —3.5t. (a) pop- 
ulation transfers between the two-qubit states |01) and |10), and (b) the corresponding adiabatic parameter 
during the passages, (c) The final population of the state |10) at a defined time t = 200 ns varies with the 
dissipation 7 = FT (with T = 4 x 10^^). The red line in (c) is obtained by numerical solution to the 
Schrodinger equation related to the Hamiltonian ^ and the dotted blue line is the analytic solution to the 
dynamics for the reduced Hamiltonian Q. (d) Probability of transfer from the states 1 01) to 1 10) varies with 
the dissipation parameter 7 and the evolution time t. 

e.g., 2ms (if it is still shorter than the decoherence time of the system). However, for the strong 
dissipations, 7 ~ 1 and 7^1, the population transfer should be achieved within sufficiently-short 
time interval. 

In summary we have investigated the Stark-chirped rapid adiabatic passage (SCRAP) of a 
driven dissipative two-level system. As a simplified model, we describe the dissipation of the 
system by adding a phenomenal parameter F to the chirped excited state of the system. Then, by 
solving the relevant Schrodinger equation we then discuss how the practically-existing dissipation 
influences the population transfer between the two selected levels of the system. We have found 
that the desired SCRAP probability is related to the effective dissipative parameter 7 = FT (with 
T being the time interval of population passage), and consequently we can divided the dissipation 



into three regions; (i) weak dissipation (7 ^ 1), (ii) strong dissipation (7 ~ 1) and (iii) very strong 
dissipation (7 ^ 1). Our numerical results clearly show that, in the weak dissipation regime, the 
SCRAP-based quantum computing scheme still works well; while in the strong dissipation regime 
the fidelity of quantum gate implemented by the SCRAP technique decreases manifestly. Cer- 
tainly, if the system works in the very strong dissipation regime, then the SCRAP technique can 
not be utilized to implement quantum computing. 

Our generic discussion has been demonstrated with a typical quantum computing system, i.e., 
the flux-biased Josephson qubits. In this specific model we have found that the loss of the transfer 
efficiency of the SCRAP is related to both the non-adiabatic transitions between the adiabatic 
passage paths and the decay of the adiabatic passage paths. During the passage for implementing 
the single-qubit gate, we find that the loss owing to dissipation-induced transition between two 
adiabatic passage paths is really small and thus negligible. For the two-qubit gate, we find that the 
dissipation-induced transition between two adiabatic passage paths vanish, and only the decay of 
the adiabatic passage paths exists. Based on this analysis we have delivered a proper approach to 
implement the quantum logic gates in such a system in the presence of dissipation. Our results 
provide quantitative estimates of the population losses during the SCRAPs, and thus should be 
useful for the realistic qubit operations. 
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